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2D  WTMM  Methodology  :  PhD  work  of  N.  Decoster 
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WTMM  Methodology  :  Skeleton 


WTMM  Chains 


WTMMM 


WTMM  Chains  at  3  different  scales 


Maxima  lines  :  WT  skeleton 


Local  roughness  characterization  :  Holder  exponent 


/(xo  +  Au)  -  /(xo)  ~  A'*<^«^(/(xo  +  u)  -  /(xo)) 


•  Monofractal  Image 


•  Multifractat  Image 
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WTMM  Method  :  multifractal  formalism 


Singularity  spectrum  ; 

D{h)  =  dnir  G  R‘^,  h{r)  =  h} 


Analogy  with  statistical  physics  :  compute 

partition  functions 

Z{q,a)  =  (M^(r,o))’  ~ 

C{a) 


Legendre  transform  : 

D{h)  =  ining(^h  —  r(g)) 


log2(a) 


WTMM  Method  :  multifractal  formalism 


Singularity  spectrum  ; 

D{h)  =  dnir  G  R‘^,  h{r)  =  h} 


Legendre  transform  : 

D{h)  =  ining(^h  —  T(q)j 


Analogy  with  statistical  physics  :  compute 

partition  functions 

Z(q,a)  =  (M^(r,o))’  ~ 

C{a) 

7i{q^a)  =  ^  In  lMi/,(r,  a)|  W^{r,a)  ~ 

r(a) 

'D{q^  a)  =  ^  In  |  W^(r,  a)|  W^(r,  a)  ~ 

z:(a) 


log8{Z(a,q)) 


Application  to  synthetic  monofractal  surfaces 


fractional  Brownian  surfaces  :  (r) 

^  H  <  0.5  :  anti-correlated  increments 
^  H  =  0.5  :  uncorrelated  increments 
^  H  >  0.5  :  correlated  increments 


H=l/3 


Theoretical  Predictions  ; 

T{q)  is  linear : 

T{q)  =qH-2 

multifractal  spectrum  is 
degenerated  : 

D(h  =  H)  =  2 


Application  to  synthetic  multifractal  surfaces 


Multifractal  (Fractionally  Integrated 
Singular  Cascades)  surfaces 


FISC 


Theoretical  predictions  : 

r(<7)  is  non-linear 
T(q)  =  -2-q{\-H*) 

-  log2(Pl^  +  P2^) 
singularity  spectrum  is  a  non¬ 
degenerated  convex  curve 


Mammography  and  breast  anatomy 


Goals  :  using  WTMM  method  to  help  in  diagnosis 

of  breast  cancer 


Fkj.  a.  ^  a  nalomio 


Fig,  1-1  ij. "  tjmt§  tHictufii-ioiMtabv,  seton  We^kigs. 


What  is  breast  cancer  ? 


M  malignant  tumor  of  mammal  gland 
3  incidence  :  30000  new  case  each  year  in  France 
3  prevention  is  very  difficult  (as  opposed  to  lung  cancer) 
3  hereditarity  :  5  to  10  %  only  (BRCA1/2  genes) 

3  forecast  depends  on  the  tumoral  volume  at  diagnosis 

SCREENING  using  mammography 


Radiological  anomalies 


distortions 


Digitalized  mammographies  :  texture  analysis 


Ji  dense  breasts  :  more  difficult  to  diagnose 
3  only  2  classes  of  monofractal  properties 

Digital  Database  for  Screening  Mammography: 

http : / /marathon . csee . usf . edu/Mammography /Database . html 

Dense  breast  Fatty  breast 


log8(Z(a.q)) 


Application  of  2D  WTMM  methodology  in  mammography 


Tissue  classification  :  dense 


Dense  breast : 
monofractal,  H  —  0.65 
persitent  correlations 


Application  of  2D  WTMM  methodology  in  mammography 


Tissue  classification  :  dense  vs  fatty 


:(a) 


30  - 
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^  Dense  breast : 

monofractal,  H  =  0.65 
persitent  correlations 

^  Fatty  breast : 

monofractal,  H  =  0.30 
anti-persitent  correlations 


Application  to  digitalized  mammographies 


Colored  Maps  : 

segmentation  of  dense  h  >  0.52  areas  and  fatty  h  <  0.38  areas 


Application  to  digitalized  mammographies 


Colored  Maps  : 

segmentation  of  dense  fi  >  0.52  areas  and  fatty  h  <  0.38  areas 


Microcalcifications  detection 


Segmentation  of  WT  skeleton  lines  : 
microcalcifications  vs  background  texture 


Background  lines 
Microcalcifications 

almost-punctual  objects  behave  like 
’Dirac’  shapes  (h  =  —1) 


Cluster  of  microcalcifications 


study  of  microcalcification  spatial  distribution 

Partition  functions  : 


Djr  =  1.3 


observation  : 
fractal  ramification  of 
cluster  of 

microcalcifications 
{1  <  Dp  <  2)  seems  to  be 
correlated  to  the 
pathology’s  malignancy 


Conclusions  and  prospects  (1) 


^  the  2D  WTMM  method  provides  a  framework  for  an  automated 
measure  of  the  breast  radio-density  and  for  studying  the  fractal 
geometry  of  clusters  of  microcalcifications. 

^  further  study  is  necessary  to  validate  quantitatively  how  far 
measuring  the  fractal  dimension  Dp  could  improve 
computer-aided  diagnosis  systems  benign/malignant 


3D  scalar  WTMM  method 


3D  data 


Smoothing 

scale  a 
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Gradient  V 

- ^ 
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Wavelet  Transform 
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3D  scalar  WTMM  method  :  skeleton 


WTMM  surfaces 


WTMMM  points 


WTMM  surfaces  at  3  different  scales 


Maxima  lines  :  WT  Skeleton  (projection  along  z] 


log2(Z(q.a)) 


Test-application  to  synthetic  3D  monofractal  fields 


fractional  Brownian  fields  : 

^  H  a  0.5  :  anti-correlated  increments 
^  H  =  0,5  :  uncorrelated  increments 
»+  i?  >  0.5  :  correlated  increments 


Theoretical  predictions  : 

^  T{q)  is  linear: 

T{q)  =qH-3 

^  multifractal  spectrum  is 
degenerated: 

D{h  =  H)  =  3 


log2(a^(q,a)) 


Test-application  to  synthetic  3D  multifractal  fields 


3D  multifractal  fields  (Fractionally  Integrated 
Singular  Cascades) 


Theoretical  predictions  : 

I®’  r(q)  =  — 2  —  q{l  —  H*) 

-  log2(Pl^  +P2‘*)- 
with  Pi  +  P2  =  1 

^  singularity  spectrum  is  a 

non-degenerated  convex  curve 


3D  dissipation  field  :  isotropic  turbulence  DNS 


pseudo-spectral  code,  (512)'^  grid,  Rx  —  216  (M.  Meneguzzi) 


3D  WTMM  methodology  vs  Box-Counting  algorithms 


“Box-Counting”  algorithm,  binomial  fit  with  pi  =  0.3  and  p2  =  0.7 


Pi  P2  =  ^  •  diagnoses  a  conservative  multiplicative  structure 


dissipation 


binomial  model : 

-2  —  q  —  log2  (pi^  +P2'') 


3D  WTMM  methodology  vs  Box-Counting  algorithms 


^  “Box -Counting”  algorithm,  binomial  fit  with  pi  =  0.3  and  p2  =  0.7 

:  diagnoses  a  conservative  multiplicative  structure 


P1  +  P2  —  1 


“3D  WTMM”method  reveals  a  non-conservative  multiplicative  structure  : 
binomial  fit  with  =  0.36  and  P2  =  0.80 


Pi  +  P2  #  1 


dissipation 


binomial  model : 

T{q)  =  -2-q-log2(pi^+P2’) 


remark: 

Pi  = 

^  P1+P2 

^  inconsistent 
box-counting  ! 


Self-similar  multifractal  vector-valued  measure  (2D  case) 


Falconer  et  O’Neil  (1995) 


scalar  measure  {r  :  lim 

1^0  log  I 

Z(q,l)  =Zi^lUl) 


ol  —  h  ^ 


^  vector-valued  measure 


lim 

1^0 


logi 


a  =  h  +  2 


Self-similar  multifractal  vector-valued  measure  (2D  case) 


es-  r  (n\  = 

M  V  y  /  log  2 

'*■  D^^{h)  =  f^{cx  -  2)  =  inf^{qh  -  T^{q)). 


Tensorial  wavelet  transform  (2D  case) 

1 .  Tensorial  wavelet  transform  of  field  V  =  (X^,  V2)  : 


T^[V](b,a)  =  (T^jy,](b,a))  = 


T^AVi]  T^AV2] 
T^AVi]  T^AV2] 


3v>i[^Kb,a)  =  a  ^  ^ (r  —  b)) V5(r), j 

Direction  of  greatest  variation  of  vector  field  : 

inn  r.n,  l|Tv-[V].C|| 

|T^[V]|  =  sup - — - - 

c^o  ||C|| 

Singular  value  decomposition  of  WT  tensor: 

T^[V]  =  (g)  .( ‘"T"  <T°,„  )•  (^)^ 
Tensorial  wavelet  transform  ; 


=  1,2 


[V]  tx)  —  ^TjnaxCr^ 


Tensorial  2D  WTMM  methodology 


Data 


Tensorial  wavelet  transform 


T|/j,max[^]  (t>,  tt)  — 


max 


Modulus  Maxima  cr^ax  chains  of  tensorial  wavelet 
transform  at  scale  a  : 


(6,  a)/ 


dcr 


max 


dG 


=  0  et 


max 
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Tensorial  2D  WTMM  methodology:  Skeleton 


WTMM  Chains 


LI  V/^i  ■  JrU  !>;j 

T  cy^-  C 
■A)  J' 

P 

j.  ^  ('.U  -  C 


WTMMM 


WTMM  Chains  at  3  different  scales 


Maxima  lines  :  WT  Skeleton 


fractional  Brownian  fields  :  B//(r) 


Spectral  method  simulation 


Monofractal  2D  vector  fields 


Theoretical  predictions  : 

^  linear  T(q):  T(q)  =  qH  —  2 

^  degenerated  singularity  spectrum: 

D{h  =  H)  =2 


2D  self-similar  multifractal  vector-valued  measures 


Self-similar  multifractal  vector-valued  measures 
(Falconer  and  O’Neil’s  model) 


Theoretical  predictions: 

T{q)  =  —  log2(pi^H-p2‘^+P3‘^+P4^) 
Pi  =  P4  =  0.5,  p2  =  2  and  =  1 


m  vectorial  box-counting 
o  vectorial  2D  WTMM  method 


vectorial  box-counting  is  less  accurate 


Tensorial  3D  WTMM  method:  turbulent  velocity  field  (Rx  =  140) 


Tensorial  3D  WTMM  method:  singularity  spectrum  of  velocity 


A  ^ 

J  parabolic  fit:  x(q)  =  —  Q  —  Cig  —  C2^ 
^  intermittency  coefficient  C2  =  0.049  di  0.004 


1D  increments  method; 

longitudinal:  C2(Svl)  ^ 
transverse  :  C2(Svt)  ^ 


0.025 

0.040 


Tensorial  3D  WTMM  method:  turbulent  vorticity  field  (Rx  =  140) 


Tensorial  3D  WTMM  method:  singularity  spectrum  of  vorticity 


□  vorticity 


Tensorial  3D  WTMM  method:  singularity  spectrum  of  vorticity 


0  2  4  0  2  4  -1  -0,8  -0.6  -0.4  -0.2 

logaCa)  loge{a)  h=a-3 


□  vorticity 

o  D^{h  +  1)  spectrum  translated  velocity 
=>  same  3D  intermittency  coefficient ! 


Conclusion 


assessment: 

^  WTMM  multifractal  analysis:  moving  towards  vector  fields 
outlooks  : 

^  better  understanding  of  the  information  embedded  in  the  WT  tensor. 

^  identification  of  coherent  structures  in  turbulence  using  WT  tensor’s  smallest  singular 
value:  vorticity  fi  laments  or  sheets. 

others  applications  :  astrophysics  (interstellar  medium,  interstellar  turbulence),  MHD, 
geophysics, ... 

Thanks  : 

^  E.  Leveque,  Laboratoire  de  Physique,  ENS  Lyon  (Turbulent  flows  DNS). 
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